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The Hanbury Brown-Twiss (HBT) effect is essentially a quantum interference of one particle with 
another, as opposed to interference of a particle with itself. Conventional treatment of identical 
particles, especially when they are entangled, is marred by attaching labels to the particles which 
appears contradictory. A recently introduced label-free approach to indistinguishable particles is 
described, and is used to analyze the HBT effect. Quantum wave-packets have been used to provide 
a better understanding of the quantum interpretation of the HBT effect. The effect is demonstrated 
for two independent particles governed by Bose-Einstein or Fermi-Dirac statistics. The HBT effect 
is also analyzed for pairs of entangled particles. Surprisingly, entanglement has almost no effect on 
the interference seen in the HBT effect. In the light of the results, an old quantum optics experiment 
is reanalyzed, and it is argued that the interference seen in that experiment is not a consequence of 
non-local correlations between the photons, as is commonly believed. 

PACS numbers: 03.65.Ud, 42.50.C 


I. INTRODUCTION 

Interference of waves is a very old and well studied 
subject. Waves emanating from two sources give rise 
to interference, provided that there is a coherence in 
the phases of the two. This can be understood with 
the following argument. Suppose there are two sources 
A and B, waves emerging from which fall on a distant 
screen. At a point xi on the screen, the contribution of 
the two classical waves can be written as 

£{xi) = (1) 

where k is the wave-vector of the two waves, TAifTsi 
are the displacements as shown in Fig. 1, and (p is the 
difference in the phases of the two sources. Assuming 
n, jS to be real, the intensity at the point xi is then given 
by 

I(xi) = |£:(xi)|^ = a^+|3^+2ajScos(A:(rAi - rgi) + (p). (2) 

The above expression represents interference between 
the wave from A and B, at point xi on the screen. One can 
see that if the phase difference between the two waves, 
p, fluctuates randomly with time, which amounts to 
integrating the above over p, the interference will be 
washed away. 

(I(xi)) +jg2_ (3) 

Thus two independent, incoherent sources of waves do 
not give rise to interference. 

Hanbury Brown and Twiss carried out an experiment 
with radio waves which involved correlating intensities 
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FIG. 1. A schematic diagram for a two source interference, and 
also for a intensity correlation at two points. Sources A and B 
are stationary while the detectors at xi and X 2 are movable. 


at two points on the screen, coming from two independent 
sources [1]. Their experiment can be easily understood 
from the preceding example. An intensity correlation at 
two points xi and X 2 on the screen can be written as 

/(Xi)I(X2) = \£iXi)\^\£iX2)\^ 

= {a^ ++ 2aj3 cos{k{rAi - rsi) + </>)} 

{a^ + + 2ajS cos(/c(rA2 - rsi) + </>)} (4) 

After averaging over p, the p-dependent terms drop out, 
and one is left with 

(I(xi)I(x2)> ~ (a^ + 

+2{apf- cos(fc(r/ii - tbi - tai + J'b2))-(5) 

Surprisingly, although individual intensities do not 
show any oscillation, there is interference between the 
intensities at two different points on the screen. This 
is what is called the Hanbury Brown-Twiss (HBT) ef¬ 
fect. Random fluctuation of phases of the independent 
source has no effect on the oscillations seen in the in¬ 
tensity correlations. Notice that the visibility [2] of the 

above interference pattern is which is bounded 

from above by 1/2. 
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The HBT effect was explained using classical waves, 
but when the same experiment was proposed using light, 
there was a lot of skepticism, and people thought it 
would not work. The reason was that the HBT effect 
arose from two different waves interfering with each 
other, and many people believed that a photon inter¬ 
feres only with itself, two photons never interfere with 
each other [3] . Hanbury Brown and Twiss carried out the 
experiment with light and demonstrated the HBT effect 
[4]. The first quantum understanding of the HBT effect 
for light was given by Fano [5]. The HBT effect has now 
been demonstrated for quantum light [6]. Not only that, 
it has also been demonstrated with massive particles, 
both of bosonic [7-9] and fermionic [10] nature. 

The HBT effect for massive particles is intriguing, as 
it involves interference between two difference particles. 
Particles are seen to be bunched together without any 
interaction between them. It is a purely quantum me¬ 
chanical effect, and the classical wave explanation is only 
a part of the story. Full treatment of the HBT effect can 
of course be done using quantum many-body theory, as 
the particles are treated as being indistinguishable there. 
However, that maybe an overkill, as it is essentially a two 
particle effect. In the following we treat two massive par¬ 
ticles as traveling wave-packets, and use them to analyze 
the HBT effect. 


II. LABEL-FREE APPROACH TO IDENTICAL 
PARTICLES 

For a quantum treatment of the HBT effect, it is essen¬ 
tial to take the indistinguishability of the particles into 
account. However, treatment of identical particles in 
quantum mechanics, particularly those involving entan¬ 
glement, has been an issue which has been much debated 
[11-29]. The problem can be seen in the very basic sym¬ 
metric or anti-symmetric wavefunctions that are usually 
written, e.g. (pA(xi)(pB(x 2 )±(pB(xi)(pA(x 2 )- Although the 
particles are assumed to be indistinguishable, we put la¬ 
bels 1 and 2 on them. The problem gets compounded 
when one wants to write an entangled state for two such 
particles. Recently a new state-based approach has been 
introduced which does not label the particles [30]. In 
the following, we will explain this approach and use it 
to analyze HBT experiment with quantum particles. 

In this new approach, the basic assumption is that in 
dealing with more than one identical particles, the par¬ 
ticles cannot be individually addressed. This is in accor¬ 
dance with the quantum theory. The combined state of 
two particles, which may consist of single-particle states, 
is treated as a holistic indivisible entity [30]. One can¬ 
not ask for the form of this state. What one can ask, is 
what the two particle probabiliity amplitude of finding 
the two particles in two different states is. For example, 
if one particle is in the single-particle state |i/^) and one 
in the state |(/)), the two-particle state is represented as 
|(/), i/^). One may not talk about the "form” of this state 


in terms of single-particle states \ip} and \(p}. However, 
the probability amplitude of finding one particle in state 
\a) and the other in state |]S), i.e., in the combined state 
|a, jS), is given by 

= {a\(p){^\\p} + ri{a\ilj}{^\(p}, (6) 

where rj^ = 1. One can check that the state \(p, ip) is 
not normalized simple by replacing |n,j8) by in 

the above equation. The state |(/), i/^) has to be multiplied 
with , ^ — in order to normalize it. The probability 

amplitude of finding the particles at positions xi and X 2 , 
i.e., in the state |xi, X 2 ) is given by 

{Xl,X2\(p,lp} = {Xi\(p){X2\lp) + ri{Xi\lp}{X2\(p), (7) 

which is the familiar symmetric or antisymmetric two- 
particle wavefunction. The difference is that, in this case 
1 and 2 are not labels of particles, but corresponds to the 
two positions of a joint measurement. 

A one-particle operator A acts on the two-particle state 
in the following way: 

A\(p,ip} = \A(p,ip} + \(p,Aip}. (8) 

The expectation value of a one-particle operator, using 
(6) and (8), is given by 

(A) = ((j),ip\A\(l),ip} 

= {(p\A(j)) + {ip\Aip) 

which agrees with the expression of the conventional 
analysis. 

III. INDEPENDENT PARTICLES 

Let there be two particles described by two wave pack¬ 
ets, traveling along y-axis. The two particles emerge from 
two sources localized at positions xq and -xq. We de¬ 
scribed the two particles by two Gaussian wave packets 
of width e each, localized at xq and -xq, denoted hy\(j)A) 
and |(/)b), respectively. Since the particles are indistin¬ 
guishable, the combined initial state of the two particles 
can be written as 

\^m = \cPA,cPB) ( 10 ) 

This satisfies the essential requirement for HBT effect, 
that the particles be identical, in the quantum sense. The 
conventional two-particle wavefunction is then just the 
probability amplitude of finding one particle at Xi and 
one at X 2 , and can be written down using (6) as follows: 

lp(Xi,X2,0} = {Xi,X2\(pA,(pB} 

= + ri{Xi\(j)B){X2\(pA) 

1 / -(^1-^0)^ -(x 2 +xo)^ -(xi+xp)^ -(^2-^0)^ ^ 

yne \ j 


( 11 ) 
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where 77 = +1. The last line in the above equation 
specifies the Gaussian form of the states emerging from 

the sources A and B, namely, <PaM = exp ^ ^ arid 

(PbM = exp . For bosonic particles, the wave- 

function should be symmetric, and rj should be 1. For 
fermions, the two-particle wavefunction should be an¬ 
tisymmetric, requiring rj to be -1. We assume that the 
particles are traveling along the positive y-direction with 
a constant velocity vq. For simplicity we ignore the ex¬ 
plicit time evolution along the y-axis, and assume that 
evolution for a time f' just transports the wave packets by 
a distance I = vot\ The dispersion of the wave packets 
along the transverse x-direction is more interesting and 
may give rise to interference between wave packets. 

Notice that if one of the two sources produces a wave 
packet with an additional phase factor, say that 
phase factor will present in the both the terms in ( 11 ), 
and can be pulled out, becoming irrelevant. 

The Hamiltonian governing the time evolution is that 
of two free particles. The two-particle eigenstate will 
simply be |p, p') which means one particle has momen¬ 
tum p and the other p'. One can thus write 

(lb+ L) 

After traveling for a time f, the particles reach the screen. 
Time evolution of the initial state \ (pA/^B) can be worked 
out by introducing a complete set of two-particle mo¬ 
mentum eigenstates: 

v,v' 

= y e 2mn |p,p ) 

[{p\(pA}{p'\(pB} +ri{p\(pB}{p'\(pA)] 

= J]\Prp') 

[{p\<PA(mp'\<pB(t)) +ri{p\<pB{mp'mt))] 

= (13) 

where 

{p\(pAit)) = e-^ (p\(pA}, {p'\(pB{t)} = e-^ iplcpB)- 

(14) 

From the momentum representation of {(pA) and |(/)b), 
the position representation can be evaluated. The ampli¬ 
tude of finding the particles at positions Xi and X 2 then 
works out to be 

( -(^ 1 -^ 0 )^ -(^2+^oF -(xi+xq)^ -ix2-xo)^ \ 

e e^+iA e + 7 ]B e^+iA e , 

( 15 ) 


where A = 2htlm, and a = —j====. For simplicity we 

introduce the notation pt = p{x\,X 2 ,i). Joint probabil¬ 
ity density of finding the particles at X\ and x^ is given 
by 


\ i - 2 e^(xi-xo)^ - 2 e^(x 2 +xo)^ -le^jxi+XQ)^ - 2 e^(x 2 -xo)^ 

\pt\ — - B C e^+A^ e e^+A^ +B e^+A^ 6 e^+A^ 

no^ \ 

-(xi+xp)^ -(x2-xo)^ -(xi-xp)^ -(x2+xp)^ 

~\~p6 e'^+iA B e'^+iA B e^-iA B e^-iA 

-(xi-xp)^ -(x2+xp)^ -(xi+xp)^ -(x2-xp)^ \ 

+ 7]C e^+iA B e^+iA B e^-iA B e^-iA j ^ (16) 


where = A A- tAjA. The above simplifies to 


\^t\^ 


2 -2e'^{xj+x^+2x'^) 

- -B e4+A2 cosh 

no^ 


4e^(xi - X2)xo 

+ a2 


1 + 


/ 4A{xi-X2)xo \ \ 
\ £4+^2 j 


(17) 


The RHS of Eqn. (17) represents an interference pattern 
in the joint probability of detection of the two particles, 
with respect to the distance between the two positions 
xi and X 2 on the screen. For p = 1, it constitutes the HBT 
effect for massive particles, which obey Bose statistics. 
The same result will also apply to independent photons 
with the proviso that A = AL/tt, where A is the wave¬ 
length of the photons and L is the distance traveled by 
them in time t. Inf act, the relation A = AL In may also be 
used for massive particles in which case A will represent 
the d'Broglie wavelength of a particle. It is easy to check 
that, for \xi - X 2 |,xo ^ r^i + TAir^Bi + ^B 2 / the cosine 
term in (5) reduces to cos(4xo(xi - X 2 )/A), which is the 
same as the cosine term in (17), provided that A^. 

At this stage it may be useful to understand the phys¬ 
ical meaning various terms in the joint probability dis¬ 
tribution (16). In our usual classical way of thinking, we 
imagine that there is a possibility of the particle from 
source A reaching xi and that from source B reaching 
X 2 . The solids lines in Fig. 1 represents this possibil¬ 
ity and the first term in (16) represents its probability. 
There also exists a possibility of the particle from source 
A reaching X 2 and that from source B reaching xi . The 
dashed lines in Fig. 1 represent this possibility and the 
second term in (16) represents its probability. If two in¬ 
dependent particles can interfere with each other, there 
is also exists possibility of particles from A and B par¬ 
tially contributing to the particles detected at xi, and at 
X 2 . Since the particles are identical, we have no way of 
knowing which source they came from. In Fig. 1, the 
solid orange line and the dotted blue line represent the 
possibility of A and B contributing to the particle reach¬ 
ing xi, and the solid blue line and the dotted orange line 
represent the possibility of A and B contributing to the 
particle reaching X 2 . The last two terms in (16) represent 
this possibility. One can see that if the particles are not 
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identical, the last two terms would not be there. This ar¬ 
gument is in agreement with the fact that the HBT effect 
cannot be seen for particles which are not identical. 

The effect of the last two terms can also be interpreted 
as the interference between the processes represented 
by the solid and the dashed lines. For certain values of 
xi, X 2 it might so happen that the solid lines process and 
the dashed lines process destructively interfere. In that 
case, there will be no simultaneous detection of particles 
at all. 

The visibility of this interference pattern is 
-, which is bounded from above by 1. Con- 


cosh 


4:e^ixi-X2)xQ 

e^+A^ 


fast this with the classical HBT effect described by (5), 
where the visibility cannot be greater than 1/2. This 
implies that there are certain distances between the de¬ 
tectors 1 and 2 for which the probability of detecting 
particles simultaneously is zero! The probability of de¬ 
tecting two particles very close to each other is enhanced. 
This can be interpreted as a bunching effect. Remarkably, 
the particles tend to get bunched together even when 
there is no interaction between them. This is a purely 
quantum mechanical effect and has been experimentally 
observed in photons [6] as well as massive particles [7-9] . 

For rj = -1, the relation (17) implies that the probabil¬ 
ity of simultaneously detecting two particles very close 
to each other is nearly zero! Even when the particles 
strike the screen at random, there is a certain probability 
of two of them hitting the screen very close to each other. 
In the case r] = -1, the probability of landing very close 
to each other is even smaller than this random chance. It 
appears as if the particles are repelling each other. This 
is what is called anti-bunching effect and has been ob¬ 
served for particles following Fermi-Dirac statistics, e.g., 
electrons [10]. 


IV. ENTANGLED PARTICLES 


Let us now investigate the scenario where the parti¬ 
cles are entangled. There are certain sources of photons 
which generate photons in entangled pairs. Entangle¬ 
ment manifests itself in strong quantum correlations be¬ 
tween the two particles. To our knowledge, the effect of 
entanglement on the HBT effect has not been quantified. 
Einstein, Podolsky and Rosen (EPR) first drew attention 
to a momentum entangled state of two particles [31] 

ipxi -ipx2 -i'^XQp 

dpe ^ e ^ e ^ , (IS) 

oo 

which can be written in Dirac notation as 




\^epr}= \p}i\-p}2e * dp, 


(19) 


where labels 1 and 2 refer to particle 1 and 2, respectively. 
This state is for distinguishable particles. If one were to 


write an EPR state for identical particles, in our label-free 
approach, it would be the following 

= f \p,-p)e~^dp. (20) 

J —oo 

The amplitude of finding the particles at xi and X 2 is 
given by 


/ oo 

[{x^\p){X 2 \-p) 

oo 

-ilxQp 

+ri{xi\-p}{x2\p}]e n dp, (21) 

which will be just the symmetrizd or antisymmetrized 
form of (18). 

The problem with the EPR state (18) is that it cannot 
be normalized, and also it does not describe particles 
with varying degree of entanglement. To address these 
shortcomings, we introduce a generalized EPR state for 
identical particles 


\W) = C 



-i2xQp p^ 

dq\q + p,q - p)e ^ e e 




( 22 ) 

where q + p and q - p label single-particle momentum 
eigenstates, C is a normalization constant, and cr, Q are 
certain parameters. In the limit a, Q ^ oo the state (22) 
reduces to the EPR state (20), if -2xo here is identified 
with xq in the EPR state. 

The two-particle amplitude of finding them at xi and 
X 2 is given by 


-X 2 - 2 xof(T^ g-(xi +X2F/4Q2 
-\--qe~^Xi-X2+2xoEo^ ^-{xi+X2El^CE^ _ ^23) 

The state (23) is an extended version of the generalized 
EPR state introduced earlier [32]. It is straightforward to 
show that Q and a quantify the position and momentum 
spread of the particles in the x-direction. The interest¬ 
ing thing about this state is that for 2Q = Ija = e^|^, it 
is no longer entangled and reduces exactly to the sym¬ 
metric state (11) studied in the last section, which is a 
symmetrized or anti-symmetrized product of two Gaus- 
sians centered at xq and -xq. So the entangled state is 
essentially two shifted Gaussians entangled with each 
other. The state (23) is symmetric under the interchange 
of the two particles, thus describing bosonic particles. 

The stage is now set to study HBT effect with two 
entangled particles, described by the state (23). The two 
particles travel in the y-direction for a time t before reach¬ 
ing the screen. During this time, the states evolves in 
transverse x-direction too. As done in the last section, 
we ignore the time evolution in the y-direction, and only 
consider the evolution in the x-direction. If one is deal¬ 
ing with photons, one can use an alternative wave-packet 


nxirX2) = 


7tQ 
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evolution [33]. The state of the two particles, on reaching 
the screen (or detectors), is given by 


W{xi,X 2 rt) = Qe 

+ 7 ] exp 


ixi+X2)^ 

4Q2+)6 


(xi - X 2 - 2 xo)^ 


lexp 

l/a 2 + id 


(xi -X 2 + 2xof 


l/(j 2 + id 


(24) 


where Ci 




-1/4 


,6 = 


4ifit Im = lALju and L is the distance in the y-direction, 
traveled by the particles during time t. 

The probability density of joint detection of particles 
at xi and X 2 can now be calculated, and is given by 


-8nHxi+X2f -2iixi-X2)^+2xl)/a^ 


16Q4+62 0 1 / o ^+52 

8(xi -X2)xo/a^ 


cosh 

/ 

1 + 7]- 


1 /(j^ + 6 ^ 


cos 


/ 86 (xi-X 2 ho ^ 
\ l/cr4+62 


\ 


cosh 


8 (xi-X 2 ho/g^ 

1 /o^+d'^ 


(25) 


J/ 


The above expression closely resembles (17) derived for 
particles which are not entangled. Let us explore it in the 
situation when the entanglement between the two parti¬ 
cles is strong. Entanglement is strong when a is large and 
one can safely assume l/cr^ ^ 6 ^. In this limit, the cosine 
term becomes cos[ 8 (xi - X 2 )xo/ 6 ] = cos[4(xi - X 2 )xo/A], 
which means that the interference fringe width is the 
same as in the case of independent particles and also in 
the classical case. The Gaussian terms in (25) assume the 

-2{xi+X2)^ -2{{xi-X2)^+2x^) 

approximate form ^ 402 + 52 / 402 ^ ^ 2^2 These Gaus- 

sians represent broad profiles both in xi + X 2 and X 1 -X 2 . 
Therefore it appears that entanglement does not cause 
any suppression of HBT effect, and remains almost as it 
is for independent particles. 

In the opposite limit, i.e., when the entanglement goes 
to zero, one can write 4Q^ = Ijo'^ = le^ (say). In this 
limit, (25) exactly reduces to (17), as expected. 


V. THE GHOSH-MANDEL EXPERIMENT 

In the light of the preceding analysis, we now take 
a fresh look at an old quantum optics experiment by 
Ghosh and Mandel [34]. This experiment was among 
the category of first experiments showing spatial corre¬ 
lation of photons. A UV laser beam was incident on a 
non-linear crystal resulting in the production of a pair of 
downconverted photons. Such photons are are known 
to be entangled, and show quantum correlations. The 
two photons travel in different directions at a small an¬ 
gle with respect to each other. Two mirrors are used 
to bring the two photons together, and two detectors, 
in the detection plane, detect them in coincidence (see 


Virtual 
source A 
UV Pump “ - - - 
Laser 


Mirror 


Virtual 
source B 


Mirror 



Counter 


FIG. 2. A schematic diagram of the Ghosh-Mandel experiment 
[34]. A UV pump laser falls on a non-linear crystal and pro¬ 
duces pairs of entangled photons via spontaneous parametric 
down conversion, traveling in different directions. These pho¬ 
tons are deflected by two mirrors and recombined where two 
movable detectors detect them in coincidence. 


FIG. 2 ). While a single detector saw no interference, a 
coincident count of the two detectors, as a function of 
their relative separation showed an interference pattern. 
The visibility of interference was greater than 1/2, which 
demonstrated the non-classical nature of photons. This 
experiment has also been discussed in a textbook [35], 
and the interference pattern is believed to be a result of 
non-local quantum correlation between the two photons 
[34,35]. 

One would notice the close similarity between the 
Ghosh-Mandel experiment and our model system study¬ 
ing HBT effect with entangled particles. The two detec¬ 
tors in the Ghosh-Mandel experiment, at xi and X 2 , see 
the photons reaching them after getting deflected from 
the two mirrors. In effect they see the photons as coming 
from two spatially separated virtual sources A and B (see 
FIG. 2 ). With this recognition, the setup in FIG. 2 is vir¬ 
tually the same as that in FIG. 1, and our model system 
captures the essence of the Ghosh-Mandel experiment. 
Our analysis shows interference of entangled paricles 
and the visibility is close to 1 . 

However, the surprising part is that the same result 
is obtained when we use independent bosonic particles 
which are not entangled, as shown in section III. When 
independent bosonic particles are used, one sees an in¬ 
terference in the coincidence count as a function of the 
relative position of the detectors, and the visibility is 
close to 1. But that is the result that is obtained in the 
Ghosh-Mandel experiment too. This implies that in the 
Ghosh-Mandel experiment, if the photons pairs were not 
entangled, the result would be the same as that for entan¬ 
gled photons. Thus the effect seen in the Ghosh-Mandel 
experiment is essentially the HBT effect. As the inter¬ 
ference in the HBT effect is independent of whether the 
two particles are entangled or not, the Ghosh-Mandel 
experiment is not a demonstration of non-local quantum 
correlation between photons, as many seem to believe. 
However, it may be the first unambiguous demonstration 
of interference between two photons, which is a com¬ 
pletely non-classical effect in itself, and probably would 
not have been expected by Dirac [3]. 
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VI. CONCLUSION 


To summarize, we have used wave-packets to study 
the HBT effect in quantum particles following Bose- 
Einstein and Fermi-Dirac statistics, using a recently in¬ 
troduced label-free analysis of indistinguishable parti¬ 
cles. The bunching and anti-bunching has been demon¬ 
strated through a simple analysis. We have also analyzed 
the HBT effect for pairs of particles which are entangled 
in position and momentum through an EPR like state. 
These entangled particles also show an HBT effect which 
is not different from the HBT effect in independent par¬ 
ticles, in any noticeable way. 


We have also argued that the Ghosh-Mandel experi¬ 
ment is essentially the HBT effect with entangled parti¬ 
cles. However, the interference seen in that experiment, 
is not a consequence of any non-local correlation be¬ 
tween the two photons. Exactly the same effect would 
be observed if the photons were not entangled. 
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